The spin quantum Hall effect in unconventional superconductors 
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We study the properties of the "spin quantum Hall fluid" - a novel spin phase with quantized 
spin Hall conductance that is potentially realizable in superconducting systems with unconventional 
pairing symmetry. A simple realization is provided by a d^2_y2 + id^y superconductor which we 
argue has a dimensionless spin Hall conductance equal to two. A theory of the edge states of 
the d^2_y2 + idxy superconductor is developed. The properties of the transition to a phase with 
vanishing spin Hall conductance induced by disorder are considered. We construct a description of 
this transition in terms of a supersymmetric spin chain, and use it to numerically determine universal 
properties of the transition. We discuss various possible experimental probes of this quantum Hall 
physics. 



I. INTRODUCTION 



A remarkable property of a singlet superconductor is 
the occurrence of the phenomenon of spin-charge sep- 
aration |0,0|. The superconducting condensate may be 
viewed as a collection of spinless, charge 2e Cooper pairs 
that have Bose condensed. The spin on the other hand 
is carried entirely by the fcrmionic quasiparticle excita- 
tions which do not carry definite charge. This observa- 
tion is particularly important in the context of supercon- 
ductors with non-s-wave Cooper pairing leading possi- 
bly to quasiparticle excitations at arbitrarily low ener- 
gies. The best studied case is d^'^^y^ pairing in the high- 
Tc cuprates. The resulting superconducting state has 
gapless quasiparticle excitations which dominate the low 
temperature properties. The cuprates thus provide an 
opportunity to explore the low energy properties of a gap- 
less spin-charge separated system in dimensions greater 
than one. Recent work ^^ has pointed out the possi- 
bility of realizing a novel spin phase - the "spin metal" - 
in the cuprates in the presence of disorder. This phase is 
characterized by a non- vanishing finite spin diffusion con- 
stant and spin susceptibility at zero temperature, and is 
not known to exist in insulating Heisenberg spin models, 
fn this work, we explore another novel spin phase poten- 
tially realizable in superconducting systems - the "spin 
quantum Hall fluid" . This phase is characterized by a 
quantized value of the Hall spin conductance (analogous 
to the quantized Hall charge conductance in the integer 
quantum hall effect). 

We begin by showing that such a spin quantum Hall 
fluid phase is realized by two dimensional superconduc- 
tors with d^2_y2 +idxy symmetry. The d+id state, which 
has received a fair amount of attention recently n~M , has 
been known to possess various similarities with quantum 
Hall states though the precise characterization in terms 
of spin transport has not been pointed out before. In 
particular, it has been suggested that a transition from 
the dx2-y2 to the dx2-y2 + idxy superconductor may be 



driven by external magnetic fields 
tentially realizable in the cuprates. 



and hence is po- 



Here we first calculate the bulk spin Hall conductance 
of the d + id state and show explicitly that it is quantized 
to be equal to two (in units of the dimensionless spin con- 
ductance). We then use semiclassical arguments to show 
the existence of two spin-current carrying edge states as 
required by the quantization of the bulk Hall spin conduc- 
tance. A Hamiltonian describing the propagating edge 
modes is derived. We next consider the effects of disor- 
der on the d + id state. The quantization of the spin Hall 
conductance is robust to weak impurity scattering. How- 
ever if the impurity scattering is sufficiently strong, there 
can be a phase transition to a phase with vanishing Hall 
spin conductance. The properties of this transition are 
considered next. Ignoring the quasiparticle interactions, 
this transition is argued to be described by the critical 
point of a replica non-linear sigma model theory y] with 
a topological term which describes quasiparticle localiza- 
tion in a superconductor without time reversal but with 
spin rotation invariance (class C of Ref. Q). We then 
construct a network model [Q describing this transition, 
and show that it is identical to that simulated recently by 
Kagalovsky et al. |ll). We then motivate a description 
of this transition in terms of a supersymmetric (SUSY) 
spin chain. In contrast to the SUSY spin chain which de- 
scribes the usual integer quantum Hall transition p2| , [l3| , 
this SUSY chain has only a finite number, three, of de- 
grees of freedom at each site. This enables the efficient 
use of a numerical technique - the density matrix renor- 
malization group (DMRG) - which has been successfully 
used for accurate calculations of the properties of quan- 
tum spin chains in other situations p4| . We present nu- 
merical results for a number of universal critical proper- 
ties of the transition. Some of these have been obtained 
before from the network model simulations [0| . Very re- 
cently, Gruzberg, Ludwig and Read [051 have provided a 
mapping of this transition to classical percolation and de- 
termined exact values for various critical exponents. Our 



numerical results are in excellent agreements with these 
exact values. We conclude with a general discussion of 
various experimental probes of the physics discussed in 
this paper. 
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Equivalently, we may require that the second quantized 
Hamiltonian Ti in Eq. be invariant under 



II. BULK SPIN HALL CONDUCTANCE OF THE 
D+ID SUPERCONDUCTOR 

We begin by defining the spin Hall conductance. In 
general, the spin conductance measures the spin current 
induced in the system in response to a spatially vary- 
ing Zeeman magnetic field. The spin Hall conductance 
measures the spin current in a direction transverse to the 
direction of variation of the external Zeeman field. More 
precisely, a Zeeman field B^{y) along, say the z-direction 
of spin, which depends only on, for instance, the spatial 
y-direction leads to a current j| of the z component of 
the spin along the spatial x-direction given by 



Jx 



dy 



(1) 



with crj being the spin Hall conductance. (Note that 
the analog of the "electric" field for spin transport is the 
derivative of the Zeeman field) . Just like the usual Hall 
effect, trj is zero in the presence of parity and time re- 
versal invariances. The d + id superconductor is neither 
parity nor time reversal invariant and hence can have a 
non- vanishing a^ . 

Before proceeding further, it is worthwhile to recall 
some general properties of singlet superconductors with 
no time reversal invariance. Consider a general lattice 
BCS Hamiltonian for such a superconductor: 



^-E 



tij 2_^ '^ia^ja + ^ij'^i^C-lj + ^i 



CjiQT 



(2) 



where i,j refer to the sites of some lattice. Hermiticity 
implies tij — t*j, and spin rotation invariance requires 
Ay — Aji. 

It is often useful to use an alternate representation in 
terms of a new set of d-operators defined by: 



The Hamiltonian, Eq. 0, then takes the form 

A,: 
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-t* 



d-J = ^4Hijdj 



(3) 



(4) 



d — > iayd* . 



(7) 



The advantage of going to the d representation is that 
the Hamiltonian conserves the number of d particles. 
Note that the transformation Eq. implies that the 
number of d particles is essentially the z component of 
the physical spin density: 



Sf 



d\d, - 1 



(8) 



A spin rotation about the z axis corresponds to a C/(l) 
rotation of the d operators. This U{\) is clearly present 
in the d Hamiltonian. Invariance under spin rotations 
about the x or y axes is not manifest though. 

Now consider the particular case of a dj;2_y2 + idxy su- 
perconductor p^ . In momentum space, the Hamiltonian 
is 



n^ 



koi 



ekcl^Ckc. + (Afc4^cLfc^ + H.c.^ 



(9) 



where Cfe is the band dispersion and A/j = Aq cos(20/c) — 
iAxy si'ii{29k) with tan(0fc) — ky/kx- It is sometimes use- 
ful to think in terms of a lattice version of the d + id 
superconductor. This has been formulated by Laughlin 
H. Translating to momentum space, for a square lattice, 
we have e^ ~ {cos{kx) + cos(/Cy)), A^ ~ Ao(cos(A:a;) — 
cos{ky)) — iAxy sin(fcx) sin(fc^) which has the same sym- 
metry under four-fold rotations of the lattice as the form 
written down earlier. 

The parameter A^y measures the relative strength of 
the dxy and d^^-y^ components. A^y = corresponds 
to the familiar dx^-y-^ state. In this limit, the gap func- 
tion Afc vanishes at four points of the Fermi surface and 
there are gapless quasiparticle excitations at these four 
nodes. A low energy theory of the dx2-y2 superconduc- 
tor can be obtained H by linearizing the dispersion re- 
lation of these quasiparticles around the nodes. We put 
Ti(A:) = Cfc, T2{k) = i(7yC_j, for ky > to write 



n = Y^T\k){ekT, + AkTx)T{k) 



(10) 



Writing tij 

real symmetric and 6^ 



afj+ibij, Ay 



ia^! with Si 



^ji-! 



'bji, real antisymmetric gives 



J^ij — ^f^ij 



a-ij ■ cr, 



(5) 



where Si are the three Pauli matrices. Note that SU{2) 
spin rotational invariance requires 



where the prime indicates a sum over ky > and r are 
the Pauli matrices in Ti,T2 (particle- hole) space. If 
{Ki,K2) are the two nodal directions with ky > 0, we 
may just keep modes near {Ki,K2). Linearizing ek and 
Afc near the nodes, we get the following low energy theory 
for the dx2-y2 superconductor: 



rL= f (fXtpli-ivpdxT^ 



(11) 



Here X = -^{x + y) and Y = -yK{-x + y). The field i/i^ is 



V2 



V2' 



the Fourier transform of 'ipi{k) = T{Ki + k) for i = 1,2. 
Each ipi thus has four components if^iaa where a is the 
particle-hole index and a the spin index. The V'l trans- 
form as spinors under SU{2) spin rotations. This Hamil- 
tonian is manifestly invariant under spin SU{2). (It also 
has additional U{1) symmetries that can be related to 
momentum conservation that holds in clean systems [p[). 
The physical charge density is of course not conserved as 
is already apparent from Eq. 0. 

It is useful at this stage to express the original real 
space electron operators in terms of the low energy con- 
tinuum fields. This is easily seen to be 
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'jiT 
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AKa -X 



V'jU + e" 



iKa -x 



^j2V 



(12) 



(13) 



with a sum over the node index (j = 1, 2) understood. 

Now consider introducing a small idxy component, i.e 
letting Aa;j, > 0. For small l^xy^ we may work with the 
low energy theory Eq. O near the nodes of the d-SC. 
The id perturbation adds to the low energy Hamiltonian 
Eq. O the following term 

Um = I d^XAxy{-4>lTy^i - ^\Ty^2) (14) 

Note that this is basically a mass term for the two Dirac 
theories describing the two nodes. 

The spin density can be expressed in terms of the con- 
tinuum fields as 



S ^ -ip^ail) 



(15) 



Similarly the spin currents may also be obtained from 
Noether's theorem. 

We now perform the continuum version of the trans- 
formation Eq. ^ by defining new fields Xiaa through 



V'ioT = Xia-\ , 



(16) 
(17) 



The form of the Hamiltonian Eqs. (O) and (14) is un- 
changed under the transformation to the x fields. It is 
clear that the z-component of the physical spin density is 
essentially the density of the x particles. A spin rotation 
about the z axis corresponds to a U{1) rotation of the x 
fields. This U{1) is clearly present in the x Hamiltonian. 
Once again invariance under spin rotations about the x 
or y axes is not manifest. 



The d operator in real space may also be expressed in 
terms of these continuum fields as 



c^{x) = d^{x)^e''^J-'^Xji\ 



-iKi -x 



X]2[, 



c\{x) = di{x) 



-iK-j -x 



XjH 



e'^^-^X,2T, 



(18) 



(19) 



with a sum over the node index {j = 1,2) understood. 
Note that the symmetry transformation Eq. M implies 
symmetry of the Hamiltonian under 



Xjaa -> i{<^y)a0x]al3- 



(20) 



The calculation of the spin Hall conductance is simpli- 
fied by choosing the external Zeeman field to be oriented 
along the z spin direction. In that case, the spin Hall 
conductance is just the charge Hall conductance of the x 
fields. The result is well-known [|l8[: The contribution of 
each Dirac species is 

2Sgn(A.,)i^ 

We have introduced the quantum of spin conductance 
^ J^^ = i^- ^^ there are now four Dirac species, we 
obtain for the spin Hall conductance (in units of ^) of 
the d + id superconductor: 



a' = 2 sgniAxy). 



(21) 



This is the main result of this section. (If we repeat 
the calculation for a dx2^y2 +is superconductor, we find 
a^ = consistent with the analysis in the following sec- 
tion on edge states). 

The explicit calculation above was restricted to 
\Axy\ ^ Aq. However the result Eq. |2l| holds even 
away from this limit. This is because the system is in 
the same phase for any finite non-zero value of the ratio 
-^. The quantized value of the spin Hall conductance 
is a universal property of this phase. A topological in- 
variant characterizing the d + id phase has previously 
been discussed by Volovik [|l6[. The results of this sec- 
tion provided a physical interpretation of this topological 
invariance in terms of the quantization of the spin Hall 
conductance. 



III. EDGE STATES 

A. Semiclassical argument 

As is well known from the theory of the quantum Hall 
effect, the quantization of the bulk spin Hall conductance 
implies the existence, for a system with a boundary, of 
spin-current carrying states at the edge. In particular 
a^y — 2 implies the existence of two such edge modes. 
Consider the d+id superconductor with a boundary, and 
a particle incident on the boundary with wave vector fci 



directed 45 degrees to the normal. This particle is re- 
flected to a state with a wave vector k2 also at 45 degrees 
to the normal. This particle can now Andreev reflect off 
the bulk of the superconductor and return as a hole (see 
Fig. nl). The hole moves on the reverse trajectory until 
it is Andreev reflected from the bulk back as a particle 
at wavevector ki. 




FIG. 1. Semiclassical trajectory leading to a surface bound 
state 



If the direction of fci corresponds to an angle 6i, the 
direction of ^2 corresponds to angle 02 = 9i ± ^. For 
the d + id gap Aj. — Aq cos(20/j ) — iA^y sin(20fc ) . There- 
fore one has A^j — — A^^. Thus there is a relative phase 
shift of TT for Andreev reflection at fci and ^2 respectively. 
The problem is then formally identical to that of an SNS 
junction with a phase shift of n between the two super- 
conductors. It is well known that in such a system there 
exists a state at zero energy bound in the normal layer. 
A similar situation obtains if the incident particle is at 
wave vector — fc2 when again the angle of incidence is 45 
degrees. For all other angles of incidence, the phase shift 
for the two Andreev reflections is different from tt, and 
there is no bound state. Thus there are precisely two 
surface bound states for every surface orientation of the 
d+id state. This is entirely consistent with the quantiza- 
tion of the bulk spin hall conductance to be two. This is, 
however, to be contrasted with the dj.2_y2 superconduc- 
tor where the existence of such zero energy surface states 
depends sensitively on the orientation of the interface [g| . 
Note also that for a dx2-y2 + is superconductor, there is 
no orientation of the interface for which the phase shift 
for the two Andreev reflections is tt - hence there are no 
surface bound states again consistent with the absence of 
a quantized spin Hall conductance. 

This semiclassical argument can be made precise by 
solving the Bogoliubov-de Gennes (B-dG) equations for 
the d + id superconductor in the presence of a boundary 
in the Andreev approximation. We remind the reader 
that the B-dG equations are just the eigenvalue equa- 
tions for the d-particle wavefunctions. As the calcula- 
tions are straightforward, and are very similar to those 
in the literature for the dx2^y2 superconductor, we will 
not present them here. Instead, we will show how the 
edge modes may be obtained from the continuum theory 



described in the previous section. 



B. Continuum Dirac theory 

To show the existence of edge states within the effective 
low energy Dirac theory, it is necessary that the incident 
and reflected modes (at 45 degrees with respect to the 
edge) lie along directions in momentum space which pass 
close to the nodes of the dx2-y2 order parameter. If this 
is not the case, a description of the edge states requires 
retaining bulk modes at high energies of order Aq. To 
this end, we consider an edge parallel to the y axis lo- 
cated at a; = 0. It is convenient to first rewrite the Dirac 
Hamiltonian in the original spatial coordinates {x,y): 



n 



d^x x\[-ivT''d^^ + -iiv^r"" + v,T^)dy - A^yTy]xi 



+ (1 ^ 2, a; ^ —X, Axy 



(22) 



Here we have performed a rotation about the t^ axis by 
an angle 6 = avctan^VF/vA) and defined v^ = —v cos{29) 
and Vz — vsin{29) with u^ = {vp + v\)/2. 

To establish the appropriate boundary conditions on 
the X fields at x = 0, it is necessary to use Eqs. ( p^ ) 
and ( pj| ) re-expressing them in terms of the underlying 
electron fields. As emphasized in the previous section, 
re-expressing the original BCS Hamiltonian in terms of 
the d-fermions eliminates all anomalous terms, reflect- 
ing the conservation of spin, even in the presence of the 
edge. The appropriate boundary condition is thus simply 
da{x = 0,y) = 0, which corresponds to the condition: 



Xlaaix = 0,2/) = -X2aa{x = 0, J/) 



(23) 



on the Dirac fields. 

To search for a zero energy edge state it is necessary 
to solve the wave equation which follows from the Dirac 
theory: 



i-ivT^d^ - A^yry^ix) = 0, 



(24) 



where we have assumed the (two-component) wave func- 
tion (j)a{x) is independent oi y - the coordinate along 
the edge. The appropriate solution which decays into 
the sample for a; > is readily found: (t>aix) = 
Sal Gxp{—Axyx/v). At low energies below A^y, the Dirac 
fields can be expanded in terms of this wavefunction as: 



Xjaaix,y) = {-ly 



^xy 



4>a{x)Xea{y)-, 



(25) 



with a two-component edge Fermion field, Xe{y)- Here, 
the (—1)-' factor has been included to satisfy the bound- 
ary conditions on x{x = 0, j/), and the prefactor under the 
square root has been chosen so that the one-dimensional 
edge field satisfies canonical anticommutation relations. 
The effective edge Hamiltonian can be readily obtained 



by inserting this expansion into the Dirac form in Eq. \21 
After performing the x-integration one finds: 



Hedge = dy x\a{-iVedy)Xe 



(26) 



with edge velocity Ve = wsin(26'). For the isotropic case 
this imphes Ve — vp = v^. 

The edge Hamiltonian describes a two-component one- 
dimensional chiral Fermion. Each edge mode contributes 
unity to the dimensionless Hall conductance, giving 



xy 



2. Since the edge density operator, xtXe, is pro- 
portional to the z-coniponent of spin, this is actually the 
spin Hall conductance, discussed in the previous section. 
Rotational invariance of the electron spin requires that 
the Hamiltonian be invariant under, x 
alently. 



iUyX^, or equiv- 



Xe 



icTyxl- 



(27) 



The edge Hamiltonian, Tiedge, is seen to satisfy this sym- 
metry. It is instructive to rewrite the edge Hamiltonian 
back in terms of the original Dirac fields, ipa, which trans- 
form as spinors under SU(2) rotations. In terms of one- 
dimensional "edge" Dirac fields defined via. 



■0eT = Xef, Ipel = xle^ 

the edge Hamiltonian takes the same form; 

T-iedge ^ dy 'lpl{-iVedy)iJe, 



(28) 



(29) 



with an implicit sum on a. This form is clearly seen 
to be invariant under SU(2) rotations: 4'e —^ U^'e, with 
U = cxp{iO ■ cr). 

Rather surprisingly, though, the edge Hamiltonian ac- 
tually is seen to have an additional U(l) symmetry; 
ipe — > exp(z0o)'0e- This additional symmetry can be 
traced to the conserved U(l) "charge" of the Dirac ip 
particles - called nodons in Ref. Q. Physically, this U(l) 
symmetry reflects the fact that the original BCS Hamil- 
tonian conserves the difference between the number of 
electrons at one node, say at Kj, and the node with 
opposite momentum, Kj. In the presence of impu- 
rities which break momentum conservation, this addi- 
tional U(l) symmetry will not be preserved. To see this, 
consider adding scattering impurities to the above edge 
Hamiltonians. For impurities which do not break spin 
rotational invariance, the edge Hamiltonian must still be 
invariant under, Xe ^ '^'^yxti and, moreover, conserve 
the z-component of spin: xlxe- A general form satisfy- 
ing these requirements is. 



Hi 



dy xliviy) ■f^)Xe 



(30) 



where r](y) are real functions, random in the spatial co- 
ordinate along the edge. Rewritten in terms of the -ip 
fields these become: 



Hirnp ^ dy [rj'^lpeO' 



Ipe + H.C. + Tl^Tpllp, 



(31) 



with rj — rj-j; zL irjy and a = {ax ± iay)/2. Although 
still invariant under SU(2) spin rotations ^e —^ Uipe, the 
additional U(l) symmetry is clearly not present anymore, 
due to the anomalous {4>eO'~4'e and ipla'^ipl) terms. 

Although the random terms explicitly break the U(l) 
symmetry, there is still another hidden U(l) symmetry, 
which can be revealed by making a clever change of vari- 
ables. Specifically, consider defining new fields. 



Ty exp 



— / dy' ri{y') ■ a 



(32) 



where Ty denotes a "time-ordering" along the spatial 
coordinate y. This effectively gauges away the ran- 
dom terms, and the full Hamiltonian when expressed in 
terms of the new ipe fields exhibits the U(l) symmetry 
"06 — *■ exp(i0o)'0e- This SU(2) gauge transformation will 
play an important role in analyzing the network model 
studied in the next section. 



IV. DISORDER EFFECTS 



A. Phase diagram 



We now move on to consider the effects of impurities 
on the d + id superconductor. As shown in the previ- 
ous section, the edge modes are robust to weak impurity 
scattering - hence so is the quantization of the bulk spin 
Hall conductance. Strong impurity scattering can how- 
ever lead to a transition to a phase with zero Hall con- 
ductance. It is useful to consider a phase diagram of the 
system as a function of A^y and disorder D. The general 
topology of such a phase diagram is shown in Fig. 0. 




FIG. 2. Schematic phase diagram as a function of A^y and 
disorder D; SQHF refers to the Spin Quantum Hall Fluid 



At zero D, a^y = 2 sgn{Aj;y). This spin quantum hall 
phase is stable to weak disorder as seen above. The line 
Axy = is of course the dx2_y2 superconductor. Turning 
on disorder at zero A^ localizes the quasiparticle states 
at the Fermi energy leading to a spin insulator. This 
phase should be robust to turning on a small A^y This 
is particularly clear in the lattice version of the d + id su- 
perconductor in terms of the d operators. The A^y sim- 
ply corresponds to a diagonal hopping term, and hence 
is innocuous, if weak, in a localized phase. It is clear 
then that there must be two transition lines emerging 
from the D = A^y = point (symmetrically about the 
^xy = line) separating the two quantum Hall phases 
(with a^y = ±2) from the spin insulator with cr* = 0. 

Note that the jump in a^y is by two |2y] - this is prohib- 
ited in generic non-interacting models of quantum Hall 
systems but is allowed here due to the special extra SU{2) 
symmetry. All phases have zero longitudinal spin con- 
ductance. It is interesting to ask about the behavior of 
the bulk quasiparticle density of states (DOS) p{E) as 
a function of energy in various regions of the phase dia- 
gram. It is known Q that in the spin insulator without 
time reversal invariance, p{E) actually vanishes as E'^ 
at low energies. In the d + id superconductor, for weak 
disorder, standard arguments suggest the development of 
exponentially small tails in the density of states leading 
to a weak filling in of the gap. However at disorder strong 
enough to be near the transition, we expect a larger den- 
sity of states that nevertheless vanishes on approaching 
zero energy as E"^. 



P(E) 



implies the existence of two extended unoccupied states 
at positive energies Ei and E2- These states, if filled, 
contribute —1 each to a^y. Thus as we move up in en- 
ergy and pass Ei, crj jumps from 2 to 1 and finally, as 
we pass E2, from 1 to zero. As the disorder increases and 
we approach the transition, Ei and E2 collapse towards 
zero. A nice way to move up (or down) in energy is by 
turning on an external Zeeman field as this acts exactly 
like a chemical potential for the d particles. In particu- 
lar, at finite Zeeman field, the transition splits into two 
separate ones with a^ jumping by one at each. We show 
in Fig. the phase diagram in the presence of a Zeeman 
field. 
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FIG. 4. Schematic phase diagram as a function of external 
Zeeman field B (or energy E) and a parameter 5 measuring 
the distance (at zero field) from the zero-field phase boundary 
for the 0^2 transition. 



B. Delocalization Transition 



FIG. 3. Density of states of the d particles showing posi- 
tions of extended states 



A different perspective on the phase diagram is pro- 
vided by considering the properties of the wave functions 
of the single particle Hamiltonian for the d particles. In 
the spin quantum Hall phase, a^ — 2 implies the ex- 
istence of precisely two extended states below the Fermi 
energy (each contributing unity to cr'^y). These two states 
will be at two different energies, say —Ei and —£'2 (see 
Fig. a). The particle- hole symmetry of the d Hamilto- 
nian in Eq. H {i.e. the SU{2) spin rotation invariance) 



Let us now consider the properties of the transition 
(in zero Zeeman field) in some more detail. This is a 
quantum Hall plateau transition where (tJ„ jumps by 
two. This is a new universality class for a quantum 
Hall localization transition distinct from the usual one 
described (for instance) by the Chalkcr-Coddington net- 
work model. A field theoretic description of this critical 
point in two dimensional superconductors without time 
reversal but with spin rotation symmetry is obtained on 
examining the nonlinear sigma model appropriate for de- 
scribing quasiparticle localization in such a system. In 
a replica formalism, this is a sigma model on the space 
Sp{2n)/U{n) ||,|,|ll. This field theory admits a topolog- 



ical term [|| as Il2{Sp{2n) /U (n)) = Z is non-trivial. We 
expect by analogy to the reasoning for the conventional 
integer quantum Hall transition that the sigma model 
supplemented with the topological term has a critical 
point which describes the spin quantum Hall transition. 
Introducing a Zeeman field induces a crossover to the 
conventional universality class 0] . This is of course con- 
sistent with the transition splitting into two as jumps of 
a^y by more than one are prohibited in that case. There 
is however another very significant difference between the 
spin quantum Hall transition and the conventional one. 
As mentioned above, the density of states actually van- 
ishes (at zero energy) on either side of the transition. By 
continuity, we expect that the density of states vanishes 
at the critical point as well. 

We may now formulate scaling hypotheses for various 
physical quantities of interest near the transition. On ap- 
proaching the critical point (at zero Zeeman field) by tun- 
ing the disorder D, for instance, the localization length 
^ (at zero energy) diverges as 



£,^6-". 



(33) 



where 6 is the distance from the phase boundary. Moving 
away from the critical point by turning on a Zeeman field 
also introduces a finite localization length £b diverging 



^E 



B- 



(34) 



We may now obtain, for instance, the behavior of the 
density of states p{E) at the critical point. To that end, 
note that moving away from zero energy is the same per- 
turbation as turning on a Zeeman field. Consequently, 
the localization length as a function of energy diverges 
as £,E ~ E^"^^ . The density of states may now be ob- 
tained by hyperscaling: 



piE) 



1 



Ea 



E-iUB- 



For (5 7^ 0, p{E) satisfies the scaling form 



p{E,5)r^E^''''-^Y{E5 



(35) 



(36) 



The universal scaling function Y satisfies 

Y{x ^ oo) = 1, (37) 



Y{x -^ 0) 



^3-2i/B 



(38) 



where the second line follows from requiring that p{E) 
vanishes as E^ off criticality. 



the universal critical properties. If we think of the links 
of the network model as corresponding to internal edge 
states of puddles of the quantum Hall fluid immersed in 
the spin insulator phase, then it is clear that we need 
to have two channels of propagation on each link. The 
link amplitude is the amplitude of propagation of the two 
channels. As the Hamiltonian H describing the dynam- 
ics of the system has the symmetry ayH*ay = —H, it is 
clear that the unitary time evolution operator U = e^'^' 
satisfies U'^ayU — ay. Upon restriction to a subspace 
with 2N states, this unitary operator can be represented 
by a matrix belonging to the group Sp{2N) (which is de- 
fined precisely as a 2A'' x 27V unitary matrix satisfying 
U^ayll = ay). Thus for the case of two channels, the 
amplitude for propagation is a 2 x 2 matrix belonging 
to the group Sp{2) = SU{2). The other ingredient in 
the network model is the matrix at the node connecting 
four links. Formally, this is a scattering event with four 
incoming channels and four outgoing channels. The cor- 
responding scattering matrix thus belongs to the group 
Sp{4:). Taking the link and node scattering matrices to 
be random and belonging to Sp{2) and Sp{A) respectively 
then completes the specification of the network model. 

In some recent work, Kagalovsky et al. lilt] have sim- 
ulated a network model with these symmetries and ob- 
tained numerical estimates of various critical exponents. 
Here however we will follow a different route. We will 
use the network model to motivate the construction of a 
supersymmetric quantum spin chain which can be used 
to calculate various disorder averaged properties of the 
system. For that purpose, it is actually more useful to 
consider an anisotropic version of the network model in 
which we view it as a collection of counter-propagating 
edge modes along the y-direction. Two adjacent modes 
are connected by random tunneling. (An alternative ap- 
proach to deriving a superspin chain is discussed in Ref. 
|15|.) As shown in the previous section, each edge mode 
is described by a two component chiral fermion and is 
described by the Hamiltonian 



(-ly / dy xUy) Hdy + ri,{y) ■ a) xM- 



(39) 



Here Xj refers to the j'th edge mode. The ^j(y) repre- 
sent the randomness on the links of the network model. 
To complete this Hamiltonian description of the network 
model, we need to introduce random tunneling between 
neighboring counter-propagating edge modes. The most 
general term consistent with the symmetries required of 
the Hamiltonian are 



C. Network model 

Just as for the conventional quantum Hall transition, 
it is possible to construct a network model to describe 



Y^ / dy \^~it°j{y)[x]+i{y)x3{y) ~ x]{y)X]+i{y)] 

+tj (y) ■ [x]+i {y)$X3 (y) + x]'^Xj+i]} m 



Here t^{y) and tj{y) are random variables with zero 
mean. 

Precisely this Hamiltonian for the case of just two 
neighboring edge modes has been studied in detail in 
Ref . M . It was shown that averages of physical quantities 
like the density of states and diffusion propagator could 
be obtained from an equivalent supersymmetric quantum 
mechanical problem defined by the nonhermitian "Hamil- 
tonian" : 



hff = J \{fl<^yflKA<^yA) 



mlp. 



2J2 



+ (/ia,/i)(/2a,/2) + 2(/T/i 

hbb=2Jib\bi + l)iblb2 + l) 

hfb - 2J [{b\ayfl){f^aybl) + {b,ayh){h<7yb2) 
~{flbi){flb2)~{b\h){blh)) 

K^LoUlh + b\b^ + !lh + blb2) 



(41) 

(42) 
(43) 

(44) 
(45) 



Here fj (bj) are two component fermionic (bosonic) oper- 
ators, and index j = 1, 2 labels the two edge modes. Pa- 
rameter Lo is the imaginary frequency at which we wish to 
compute averages. The constant J > is determined by 
the strength of the disorder. Its actual value is unimpor- 
tant for calculation of universal properties. We refer the 
reader to Ref. Q for further details. In the following we 
set J = 1 for convenience. This superHamiltonian gener- 
ates time evolution in the y-direction. Clearly the super- 
Hamiltonian describing the full network can be built up 
from this two edge Hamiltonian. Just as in the case of the 
superspin chain which describes the conventional quan- 
tum Hall transition |12] , the two distinct phases on either 
side of the transition correspond to the two possible ways 
of dimerizing the chain. The critical point corresponds to 
the uniform chain where the bond strength J is the same 
for all bonds. An important feature of this Hamiltonian 
that is not shared by superspin chains constructed for the 
conventional quantum Hall transition ||lj,|l^ is that the 
low energy sector of this theory is described by a finite 
on-site Hilbert space [Q of dimension D = 3: 

ID = 10) 



|2) ^ — e^pbifllO) 
|3) EE 1 e„^/t/>) . 



(46) 



This crucial simplification permits considerable numeri- 
cal and analytical progress. 

For a chain of L (even) sites, the superHamiltonian 
may be written, following the notation of Ref. O], as: 



L-2 



4 8 ■ 

^ ga s'; s';^, + i-iy V g^ s'; s^+. 



.a=l 



a— 5 



L-1 

3=0 



(47) 



Here Jj = [1 + (—1)^(5] where the relevant dimerization 
parameter J = at the critical point. We have intro- 
duced different imaginary frequencies cuj at each site to 
permit the extraction of critical properties (see the fol- 
lowing section). The constants ga are defined to be: 



' 2; a= 1,7,8 
9a = I 1; a = 3,4 

^ -2; a-2,5,6. 
In Eq. B^ we have introduced 8 spin operators: 
S' = biba + 1 S' = e^pbifl 

S = /a/a - 1 S^ = tapbafp 

S' ^ e.pflfl S' EE biU 



(48) 



(49) 



Bosonic- valued operators S*^, . . . , S''* make up the sym- 
metric sector of the Hamiltonian while fermion-valued 
operators S^ , . . . , S^ are in the antisymmetric sector. De- 
spite the fact that H is non-Hermitian, it only has real- 
valued eigenvalues. H is also defective (the left eigen- 
states or the right eigenstates do not separately span the 
whole Hilbert space), complicating the numerical prob- 
lem of diagonalizing it. 

The Hamiltonian commutes with two (fermion-valued) 
supersymmetry generators, [H, Qi] = [H, Q2] = 0, 
where 



Qi-E 



&l/.« - {-^yflb. 



ja'^ja 



Q2^Y. i-^ybU^o.+fib, 



(50) 



It is not difficult to see that the supersymmetric Hamilto- 
nian must have a unique, zero-energy, ground state. The 
right and left (ground) eigenstates are therefore annihi- 
lated by the Hamiltonian: i/|*o) = (*o|-f^ = 0. Also, 
the ground state is annihilated by the SUSY charges: 
Qil^o) — Q2|^o) — 0. All excited states appear in quar- 
tets or larger multiples of 4, half with odd total fermion 
content, and these cancel out in the partition function by 
virtue of the supertrace: 



Z = STre-^^ = Tr(-l)^/e-^^ = 1, 
where Nf is the total number of fermions. 



(51) 



V. DMRG RESULTS 

We employ the relatively simple "infinite-size" DMRG 
algorithm [H to numerically access the properties of the 
critical point 5 = 0. The fact that the ground state 
energy is exactly zero provides a valuable check on the 
accuracy of the DMRG algorithm which incurs errors 
when, as the chain length increases, the Hilbert spaces 
of the blocks grow beyond the finite limit of M states. 
Increasing Af up to limits set by machine memory and 
speed yields systematic improvement in the accuracy of 
the DMRG algorithm. In results reported below we have 
checked that M is sufficiently large to ensure adequate 
accuracy; for M > 243 there is no truncation until the 
chain exceeds length L = 12. Reasonable accuracy is 
maintained, for the case M = 256, out to L = 26: the 
ground state, when targeted, has an energy Eq which in- 
creases from zero to just Eq = 2.3 x 10~^ at L = 26. 
Furthermore, for M = 512, the ground state energy is 
only Eq = 3.2 x 10~^ at chain length L = 30, showing 
the systematic improvement in accuracy with increasing 
M. 

Reduced density matrices for the two augmented 
blocks, each of Hilbert space size D x M, are formed 
by computing a partial trace over half the chain. For the 
left half of the chain the density matrix is chosen to have 
the following symmetric form |21| : 



Py = o 



DAI 

E 

i' = l 






(52) 



a similar formula holds for the right half of the chain. 
Here vj/^, = (i,i'|*) and *fj, = (*|i,i') are, respectively, 
the real- valued matrix elements of the targeted right and 
left eigenstates projected onto a basis of states labeled 
by unprimed Roman index i which covers the left half 
of the chain and primed index i' which covers the right 
half. To compute ground state properties, '^ is selected 
to be the ground state; conversely, to find the gap, ^ is 
chosen to be one of the lowest-lying excited states. All of 
the eigenvalues of p are real and positive; these are inter- 
preted as probabilities and the {D — l)M least probable 
states are thrown away. 

To extract critical behavior, we monitor the induced 
dimerization and spin moments near the center of the 
chain as the chain length L is enlarged via the DMRG al- 
gorithm pl| ] . Dimerization is induced by the open bound- 
ary conditions as shown in Fig. 0. Spin moments are 
formed in the interior of the chain in two different ways. 
In the bulk case luj is set equal to a small, but non-zero, 
constant w > on each site, inducing non-zero spin mo- 
ments. Alternatively, the spins at the chain ends can be 
fixed by setting ujj = except at the chain ends where 
LUj is assigned a large value which completely polarizes 
the end spins, see Fig. 0. Power-law scaling of the in- 
duced dimerization and spin moments in the interior of 



the chain is expected ||2^ at the critical point 6 = 0. As 
discussed earlier, we may move off criticality either by 
dimerizing the spin chain or by turning on a finite Zee- 
man field (which is equivalent to going away from zero 
energy). There are two independent exponents related to 
these two perturbations of the critical spin chain. As in 



Section [V, we may write down scaling forms for various 



physical quantities. For a finite system size, these scaling 
forms will involve two scaling variables: the ratio -r- of 

the two localization lengths and the ratio ^. Consider, 
for instance, the density of states. This is determined by 
the boson occupancy according to M 



27: p{E) = - Re J2a + b%) (53) 



where we calculate expectation values setting loj = lu = 
—i{E + irj). Thus p{E) can be obtained from the be- 
havior of the spin operator S^ . This scales at the center 
of the chain as a function of the chain length L and the 
uniform, "bulk," imaginary frequency loj = lu as follows: 



{Sl/,)=cj^f{Lu;''-,^u;'^-) 



where the exponent 



a = 2vb — 1 



(54) 



(55) 



as required by hyperscaling (see Eq. pSJ ). When the ap- 
plied dimerization 5 = 0, this reduces to 



{Sl/^)=u;" g{Luo^-) 



L^uj 



as uo 



. 



(56) 



Here the scaling function g{x) is given, for \x\ ^ 1, by: 
g{x) = a;-"/'^^ [cix^/"^ + C2X^'''^ +...). (57) 

This equation expresses the fact that when the sys- 
tem length is much smaller than the correlation length 
{\x\ ^ 1), the DOS is an analytic, linear, function of the 
imaginary energy uj. With this scaling form we obtain 



y = 2{l- vb)Ivb 



(58) 



In what follows, we first describe the calculation of the 
exponents v., vb for the two diverging localization lengths. 
These can then be used to extract the other critical ex- 
ponents a, y using the above scaling arguments. We will 
however provide independent support for the validity of 
these scaling arguments by direct calculation. 
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Spin Moments Induced by Field Applied to Ends 



FIG. 5. Extraction of critical behavior from finite-size ef- 
fects, illustrated for the case of an ordinary quantum antifer- 
romagnetic spin chain. Dimerization of the nearest-neighbor 
spin-spin correlation function, indicated here by alternating 
strong (solid) and weak (dashed) bonds, is induced by the 
open boundary conditions. Spin moments are induced by the 
application of a magnetic field of strength ±fe to the two spins 
at the ends of the chain (shown) or formed by the application 
of a staggered field throughout the chain. 



1. Localization length exponent ly^- The localization 
length scales, as a function of the imaginary frequency to, 
with exponent vb- 



L 



(59) 



One way to determine vb is to find the crossover, for 
uniform luj = lu > 0, from power law decay of the in- 
duced dimerization to exponential decay. The induced 
dimerization at the center of the chain is defined as 



A{L,Uj) = \{SlI2-iSl/2 - Si^/2^L/2+l) 



(60) 



ea/3/i4 ^^ °^^ °^ ^^'^ ^ ^USY 



where we recall that S^ 

spin operators (each of the 7 other spin operators scale 

similarly.) It has the following asymptotic behavior: 

r C L-""; w = 
A{L,u)=l (61) 

[ C" e^"^")'^; to(w)L> 1. 

Fits to the second line in Eq. pi] permit the extraction 
of the mass gap m{uj)] then vb is determined by a power 
law fit to ■m{uj) ^ uj'^^ . We find vb = 0.55±0.1 for calcu- 
lations with M = 128, fitting over the range 20 < L < 24 
and 0.1 < a; < 0.8. A direct calculation of the gap in the 
excitation spectrum as a function of lo also yields results 
consistent with this value for vb- Note that the exact 
result of Ref. H is i/s = 4/7 = 0.5714.... 



FIG. 6. Power-law scaling of the induced dimerization, the 
bulk occupancy, and the induced occupancy with chain length 
L. In the case of the induced dimerization and the bulk oc- 
cupancy, LOj = 10~^ throughout the chain, small enough for 
the bulk occupancy to be well described by the second line 
of Eq. pq. The bulk occupancies have been multiplied by a 
factor of 10"'. The induced occupancy is obtained by setting 
ujj — everywhere except at the chain ends where it is made 



large, in this case wo = u)l- 
each of the three data sets. 



1 = 10. Straight lines are fit to 



2. Dimerization exponent v. For small luj = 10^^, fit- 
ting the induced dimerization shown in Fig. to the first 
line in Eq. ^yields a; = 1.24 ± 0.01. The dimerization 
exponent i' is related to the scaling dimension x by 



1 



(62) 



and thus i^ = 1.32 ± 0.02, close to the percolation value 
of 4/3 reported in Ref. [1^]. For the spin-1/2 Hciscn- 
berg antiferromagnet and the spin-1 antiferromagnet at 
the critical point accuracy at the few percent level was 
also achieved ||2j]. The network model simulations |ll| ] 
find 1/ ~ 1.12. Though this is close to the value we find 
numerically, and to the exact result |15|, the reason for 
the lack of more precise agreement is unclear to us. 

3. DOS exponents a and y. Drawing upon the data 
shown in Fig. |g we obtain y = 1.43 ± 0.05 by direct fit 
of the bulk occupancy at one of the central sites to the 
second line of Eq. [S^. The error is estimated by compar- 
ing results from DMRG calculations with AI = 256 and 
M = 512 and also by making power-law fits over different 
ranges of chain lengths L. This calculation of y can now 
be used to calculate i^b = 0.58 ± 0.01 in good agreement 
with the value obtained in item 1 above. 

As mentioned above, the scaling of the DOS can be 
extracted in another way: set tUj = everywhere along 



10 



the chain except at the two sites at the ends of the chain 
where it is made large. Consequently at the chain ends 
i^o) = i^L-i) = 1 but in the interior the expectation 
value {S\i^, which we call the induced occupancy, de- 
creases as the chain grows in length: 



{Si 



L/2I 



(63) 



From Fig. g we find w ~ 0.26 ± 0.02. Now, scaling re- 
lates L ^ Lu^'^'^ and hence a — wvb- Using the relation 
a = 2vb — 1, we get vb = 0.57 ± 0.02 again in agreement 
with the estimates above, and the exact result [|5|. Note 
that the density of states exponent a = 0.14 ± 0.04. 



VI. DISCUSSION 



The properties of this chiral spin metal will be quite sim- 
ilar to those of the chiral metal discussed in multilayer 
quantum Hall systems ]25| , p6[ . 

Throughout this paper, we have analyzed only the case 
of spin singlet pairing. For triplet pairing, such as in a p- 
wave superconductor, neither the spin nor the charge of 
the quasiparticles is conserved. Thermal transport still 
remains a useful way of probing quasiparticle transport. 
Arguments very similar to those used in this paper show 
that a two dimensional superconductor with Px+ipy sym- 
metry has a quantized thermal Hall conductance. For a 
layered three dimensional system, we then have a chiral 
surface sheath with diffusive thermal transport in the di- 
rection perpendicular to the layers, and ballistic thermal 
transport within each layer. Such a layered p + ip super- 
conductor is possibly realized in the material Sr2RuOi 



How may the physics discussed in this paper be probed 
if a dx2_y2 +idxy superconductor were to be found exper- 
imentally? The bulk of this paper has focused on spin 
Hall transport which is extremely difficult to measure. 
However, the thermal Hall conductance is also quantized 
in the dx^^y-i +idxy state. This can, for instance, be seen 
using the edge state theory developed in Section III . In- 



deed, if the temperature of one edge is raised by 6T rel- 
ative to the other, the excess heat current is easily seen 

27T^TSTkl 



to be 



3h 



implying a thermal Hall conductance of 



2TT^Tk% 

3h 



(64) 



-idx 



supercon- 



Thus -^^ is quantized 123] in the d^^- 
ductor. On the other hand, in the spin insulator phase, 
■^ goes to zero as the temperature goes to zero. Note 
that the charge Hall conductance is not quantized in the 
d + id phase [Q . Physically this is because any edge 
quasiparticle electrical current causes flow of supercur- 
rent in the opposite direction out to a distance of order 
the penetration depth. 

The behavior of the quasiparticle density of states may 
be probed by specific heat, spin susceptibility, or tunnel- 
ing measurements. We caution, however, that it may be 
necessary to include quasiparticle interactions, neglected 
in the theory so far, to obtain meaningful comparisons 
with experiments for these quantities. (The quantization 
of the spin and thermal Hall conductances is expected to 
be robust to inclusion of quasiparticle interactions). 

It is interesting to ask about experimental realizations 
oi d + id pairing symmetry in layered three dimensional 
superconductors. If each layer is deep in the spin quan- 
tum Hall fluid phase, then arguments similar to those for 
multilayer quantum Hall systems p^ , |26| , imply the exis- 
tence of a "chiral spin metal" phase at the surface with 
diffusive spin transport in the direction perpendicular to 
the layers and ballistic spin transport within each layer. 
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